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Abstract 

The vacuum expectation value of the surface energy-momentum tensor is evaluated for a 
scalar field obeying Robin boundary condition on a spherical brane in [D -\- l)-dimensional 
spacetime Ri x where Ri is a two-dimensional Rindler spacetime. The generalized 

zeta function technique is used in combination with the contour integral representation. 
The surface energies on separate sides of the brane contain pole and finite contributions. 
Analytic expressions for both these contributions are derived. For an inhnitely thin brane 
in odd spatial dimensions, the pole parts cancel and the total surface energy, evaluated as 
the sum of the energies on separate sides, is finite. For a minimally coupled scalar held 
the surface energy-momentum tensor corresponds to the source of the cosmological constant 
type. 
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1 Introduction 


Motivated by string/M theory, the AdS/CFT correspondence, and the hierarchy problem of 
particle physics, braneworld models were studied actively in recent years [1]. In this models, 
our universe is realized as a boundary of a higher dimensional spacetime. In particular, a well 
studied example is when the bulk is an AdS space. The problem of studying quantum effects in 
braneworld scenarios is of considerable phenomenological interest, both in particle physics and 
in cosmology. The braneworld corresponds to a manifold with dynamical boundaries and all 
fields which propagate in the bulk will give Casimir-type contributions to the vacuum energy, 
and as a result to the vacuum forces acting on the branes. In dependence of the type of 
a field and boundary conditions imposed, these forces can either stabilize or destabilize the 
braneworld. In addition, the Casimir energy gives a contribution to both the brane and bulk 
cosmological constants and, hence, has to be taken into account in the self-consistent formulation 
of the braneworld dynamics. Motivated by these, the role of quantum effects on background of 
Randall-Sundrum geometry has received a great deal of attention. The models with dS and AdS 
branes, and higher dimensional brane models are considered as well (see, for instance, references 
given in [2]). 

In view of the recent developments in braneworld scenarios, it seems interesting to gener¬ 
alize the study of quantum effects to other types of bulk spacetimes. In particular, it is of 
interest to consider non-Poincare invariant braneworlds, both to better understand the mecha¬ 
nism of localized gravity and for possible cosmological applications. Bulk geometries generated 
by higher-dimensional black holes are of special interest. In these models , the tension and the 
position of the brane are tuned in terms of black hole mass and cosmological constant and brane 
gravity trapping occurs in just the same way as in the Randall-Sundrum model. Braneworlds in 
the background of the AdS black hole were studied in [3]. Like pure AdS space the AdS black 
hole may be superstring vacuum. It is of interest to note that the phase transitions which can 
be interpreted as confinement-deconfinement transition in AdS/CFT setup may occur between 
pure AdS and AdS black hole [4]. Though, in the generic black hole background the investiga¬ 
tion of brane-induced quantum effects is technically complicated, the exact analytical results can 
be obtained in the near horizon and large mass limit when the brane is close to the black hole 
horizon. In this limit the black hole geometry may be approximated by the Rindler-like manifold 
(for some investigations of quantum effects on background of Rindler-like spacetimes see [5] and 
references therein). In the paper [6] we have investigated the Wightman function, the vacuum 
expectation values of the field square and the energy-momentum tensor for a scalar field with an 
arbitrary curvature coupling parameter for the spherical brane on the bulk Ri x S^~^, where Ri 
is a two-dimensional Rindler spacetime. This problem is also of separate interest as an example 
with gravitational and boundary-induced polarizations of the vacuum, where all calculations can 
be performed in a closed form. Note that the corresponding quantities induced by a single and 
two parallel flat branes in the bulk geometry Ri x R^~^ for both scalar and electromagnetic 
fields are investigated in [7]. For scalar fields with general curvature coupling, in Ref. [8] it has 
been shown that in the discussion of the relation between the mode sum energy, evaluated as the 
sum of the zero-point energies for each normal mode of frequency, and the volume integral of the 
renormalized energy density for the Robin parallel plates geometry it is necessary to include in 
the energy a surface term concentrated on the boundary (see also the discussion in Ref. [9]). An 
expression for the surface energy-momentum tensor for a scalar field with a general curvature 
coupling parameter in the general case of bulk and boundary geometries is derived in Ref. [10]. 
The vacuum expectation values of the surface energy-momentum tensor on the branes in AdS 
bulk are investigated in [11]. In particular, it has been shown that the surface densities induced 
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by quantum fluctuations of bulk fields can serve as a natural mechanism for the generation of 
cosmological constant in braneworld models of the Randall-Sundrum type with the value in good 
agreement with recent cosmological observations. The purpose of the present paper is to study 
the vacuum expectation value of the surface energy-momentum tensor for a scalar field obeying 
Robin boundary condition on a spherical brane on the bulk Ri x S^~^. The paper is organized 
as follows. In section 2 we consider the surface energy-momentum tensor and the eigenfunctions 
for the problem. The vacuum expectation value of the surface energy-momentum tensor in the 
R-region (the definitions of the R- and L-regions see below) are investigated in section 3. The 
corresponding quantities for the L-region are discussed in section 4. Section 5 summarizes the 
main results of the paper. 

2 Surface energy-momentum tensor 

Consider a real scalar field (p{x) on background of {D + l)-dimensional spacetime Ri x 
where Ri is a two-dimensional Rindler spacetime. The corresponding line element has the form 

( 1 ) 

with the Rindler-like (r, part and dTj'jj_i is the line element for the space with positive constant 
curvature with the Ricci scalar R = (D — 2){D — \)/r\. Line element (1) describes the near 
horizon geometry of {D + l)-dimensional topological black hole with coordinate ^ determining 
the distance from the horizon. For example, in the case of a (Z) -|- l)-dimensional Schwarzschild 
black hole one has r — rn = {D — 2)^^/4r/f, where r is the Schwarzschild radial coordinate and 
r = rn corresponds to the horizon. For the scalar field ‘^{x) with curvature coupling parameter 
C the dynamics is governed by the field equation 

-|- w? -|- QR^ (/9 = 0, (2) 

where V; is the covariant derivative operator associated with the corresponding metric tensor 
Qik- In the cases of minimally and conformally coupled scalars one has C = 0 and Q = {D — 1) /4Z), 
respectively. Our main interest in this paper will be the surface Casimir energy and stresses 
induced on a spherical brane located at ^ = a. We will assume that the field satisfies the Robin 
boundary condition 

(A,, -F nVz)99(x) = 0 (3) 

on the brane, where Ag is a constant, v} is the unit inward normal to the brane. This type 
of conditions is an extension of Dirichlet and Neumann boundary conditions and appears in a 
variety of situations, including the considerations of vacuum effects for a confined charged scalar 
field in external fields [12], spinor and gauge field theories, quantum gravity and supergravity 
[13]. Robin boundary conditions naturally arise for scalar and fermion bulk fields in the Randall- 
Sundrum model [14]. For boundary condition (3) the vacuum expectation value of the bulk 
energy-momentum tensor induced by a spherical brane is evaluated in Ref. [6]. In Ref. [10] it 
was argued that the energy-momentum tensor for a scalar field on manifolds with boundaries in 
addition to the bulk part contains a contribution located on the boundary. The surface part of 
the energy-momentum tensor is given by the formula [10] 

=5{x-dMs)nu (4) 

where the ’’one-sided” delta-function 5{x;dMs) locates this tensor on boundary dMg and 

Tik = - (2C - ll2)hik‘fn^Vnp. (5) 
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Here Kn^ is the extrinsic curvature tensor for the boundary and hik is the corresponding induced 
metric. 

Let {(/?Q,(x), (x)} be a complete set of positive and negative frequency solutions to the 

field equation (2), obeying boundary condition (3). Here a denotes a set of quantum numbers 
specifying the solution. By expanding the field operator over the eigenfunctions ipa{x), using the 
standard commutation rules, for the vacuum expectation value of the surface energy-momentum 
tensor one finds 

^Q|^(surf)|Q) = ^(^x-,dMs){0\Tik\0), (0|rifc|0) ='^Tik{<faix),(p*aix)}, (6) 

a 

where |0) is the amplitude for the vacuum state, and the bilinear form Tik{^p,'ip} on the right of 
the second formula is determined by the classical energy-momentum tensor (5). To evaluate the 
vacuum expectation value of the surface energy-momentum tensor we need the eigenfunctions 
ipQ.{x). In the consideration below we will use the hyper spherical angular coordinates = 

{01, $ 2 ,..., On, 4>) on with n = Z) — 2, 0 < < TT, A: = 1,..., n, and 0 < 0 < 27r. In these 

coordinates the variables are separated and the eigenfunctions can be written in the form 

Vc.{x) = Co,f{i)Y{mk-,^, (7) 

where nik = {mo = I, mi,... mn), and mi, m 2 , ■ ■ ■ mn are integers such that 

0 < mn-i < ■ ■ ■ <mi <1, -mn-i < mn < mn-i, (8) 

Y{mk','d,4>) is the surface harmonic of degree I [15]. The equation for /(^) is obtained from 
field equation (2). The corresponding linearly independent solutions are the Bessel modified 
functions I±iuj{^i^) and Ki^{\i^) with the imaginary order, where 

A; = —\ll{l + n) + (n{n + 1) + m‘^rjj. (9) 

rn ^ 

The eigenfrequencies are determined from the boundary condition imposed on the field at ^ = a. 
The brane divides the spacetime into two regions with ^ > a (R-region) and 0 < ^ < a (L-region). 
The vacuum properties in these regions are different and we consider them separately. 

3 Surface energy in the R-region 

For the R-region the unit normal to the boundary and nonzero components of the extrinsic 
curvature tensor have the form 

n'- = S{, Koo = a, (10) 

and /(^) = Kii^{Xi^). For a given A;a, the corresponding eigenfrequencies uj = uoj = ujj{Xia), 
j = 1, 2,..., are determined from boundary condition (3) and are solutions to the equation 

AKi^{x) + xK[^{x) = 0, x = Xia, A = Aga, (11) 

where the prime denotes the differentiation with respect to the argument of the function. For 
yls > 0 this equation has purely imaginary solutions with respect to u). To avoid the vacuum 
instability, below we will assume that Ag < 0. Under this condition all solutions to (11) are real. 
The coefficient in Eq. (7) is determined by the normalization condition. Using the relation 

j\Y{mk;^,4>)\^ d^l = N{mk) (12) 
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for the spherical harmonics (the explicit form for N(mk) will not be necessary in the following 
consideration), one finds 




~.n+l 
' H 


1 

^{ruk) -§;^Ki^{\ia)\u,=u,. 


(13) 


where for a given function F{x) we use the notation 


F{x) = AF{x) + xF'{x). 


(14) 


Substituting the eigenfunctions into the mode-sum formula (6) and using the relations 
Kiujj{Xia)Iiujj{Xia) = 1 and 


i Af(mfc) Sd' 


(15) 


the vacuum expectation value of the surface energy-momentum tensor can be presented in the 
form 


{0\rt\0) 


2r^~^aSD 


2CSf6^ + (4C - l)A5^ 


l,k = 0,2,... , D, 


(16) 


and (0|rj'^|0) = 0, with So = 27r^/^/r(Zl/2) being the total area of the surface of the unit sphere 
in H-dimensional space, and 


/r(^)=^D,E 


KiiOj (Ajfl) 


—K 

1=0 j=i 


, (A;0) 


(17) 


Here and below the quantities for the R- and L-regions are denoted by the indices R and L, 
respectively, and we use the notation 


Di = {2l + D- 2) 


T{l + D- 2) 
T{D-l)l\ 


(18) 


for the degeneracy factor. The vacuum expectation value of the surface energy-momentum 
tensor (16) has a diagonal structure: 


(Olr^'^lO) = diag (e^^\0, -p^^\ ..., , (19) 

with the surface energy density the stress 


and with the equation of state 


p(R) 


AIr{A) 

2rH~^a 


(1-4C), 


(R) _ _ 


e'= 


1 + 


2C 


yl(4C-l)J 


P 


(R) 


( 20 ) 


( 21 ) 


For a minimally coupled scalar field, the latter corresponds to a cosmological constant induced 
on the brane. Note that the vacuum expectation values of the field square on the brane is also 
expressed in terms of the function Ir{A): 


(0|7^'|0)^=a 


Ir{A) 


( 22 ) 
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The quantity (17) and, hence, the surface energy-momentum tensor diverges and needs some 
regularization. Many regularization techniques are available nowadays and, depending on the 
specific physical problem under consideration, one of them may be more suitable than the others. 
Here we will use the method which is an analog of the generalized zeta function approach. We 
define the function 

OO 

= X! ACrCS) ^la), (23) 

1=0 


where 


Cr(s, Aio) 





(24) 


Note that for Dirichlet boundary condition this function vanishes. The computation of vacuum 
expectation value for the surface energy-momentum tensor requires an analytical continuation 
of the function Tr(s) to the value s = 0, 


Ir{^) = Fr{s)\s=o- 


(25) 


The starting point of our consideration is the representation of the function (24) in terms of 
contour integral 


CR(s,a;) 


2m Jc KM 


(26) 


where C is a closed counterclockwise contour in the complex z plane enclosing all zeros ojj{x). 
The location of these zeros enables one to deform the contour C into a segment of the imaginary 
axis {—iR, iR) and a semicircle of radius R in the right half-plane. We will also assume that the 
origin is avoided by the semicircle Cp with small radius p. For sufficiently large s the integral 
over the large semicircle in (26) tends to zero in the limit R ^ oo, and the expression on the 
right can be transformed to 


Cr('S,x) 


J_ [ dzz-^^^ 

2m Jc, KM 


1 TTS 

— COS — 
TT 2 



dzz ^ 


Kz{x) 

Kz{x)' 


(27) 


Below we will consider the limit p ^ 0. In this limit the first integral vanishes in the case 
s = 0, and in the following we will concentrate on the contribution of the second integral. For 
the analytic continuation of this integral we employ the uniform asymptotic expansion of the 
MacDonald function and its derivative for large values of the order [16]. We will rewrite this 
expansion in the form 


Kz{x) ~ 



^-zri{xlz) 

(x^ -i- z^yi'^ 


OO 


E 

< 7=0 


(-l)%g(f) 

(x^ -I- ’ 


(28) 


where 


t = 


r/(x) = VT + x^ -|- In ■ 


M) = 


Uqit) 


(29) 


and the expressions for the functions Uq{t) are given in [16]. From these expressions it follows 
that the coefficients Ug{t) have the structure 


^qit) = E 

m=0 


(30) 
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with numerical coefficients Uqm- From Eq. (28) and the corresponding expansion for the deriva¬ 
tive of the MacDonald function we obtain the asymptotic expansion 


Kz{x) ~ 



-I- 

< 7=0 


(-l)^hg(f) 

{x^ -t- ’ 


(31) 


where 

= + (32) 

and the expressions for Vq{t) = X)m=o are presented in [16]. Note that the functions 

(32) have the structure 

q 

~ ^ ( Vqmt , Vqm — Vqm + A.Uq—\,m, • (33) 

m=0 

From Eqs. (28) and (31) we can find the asymptotic expansion for the ratio in the second 
integral on the right of formula (27): 


Kz{x) 

Kz{x) 


1 

(x^ -I- 


^ i-mqjt) 

^ (x^ -|- 2;2 [)ij/2 ’ 


(34) 


where the coefficients Uq{t) are defined by the relation 

-1 


E(-i) 


< 7=0 




E(-i) 


< 7=0 




= E 

< 7=0 


'fQ 


and similar to (30), (33), are polynomials in t: 


u,{t) = Y.u,it‘‘■ 

j-o 

The hrst three coefficients are given by expressions 

Uo{t) = = 




Now let us consider the function 


1 OO 

/ , 1 Tl S X—> _ 

Fa{s) = -- cos — X] A 




TT 


1=0 


KziXia) 


(35) 


(36) 


(37) 


We subtract and add to the integrand in this equation the first N terms of the corresponding 
asymptotic expansion. This allows us to split (37) into the following pieces 

Fr(s) = 4“)(s) + F«(s), (38) 
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where 


4 ”(») 


1 


- COS — y Di 



(z2 + A2a2)(<?+i)/2’ 


1 ITS 

-COS- 

TT 2 


OO 


EA 

1=0 



y(Aza) +^y„(z2 + A2a2)(</+i)/2 


(39) 

(40) 


and 

t = z/^z‘^ + \ja?. (41) 

For > H — 1 the expression for (s) is finite at s = 0 and, hence, for our aim it is sufficient 
to subtract = H — 1 asymptotic terms. At s = 0 the function (s) is finite for p = 0 and 
we can directly put this value. The integral over z in the expression for F^^\s) is finite in the 
limit p ^ 0 for 0 < Res < 1. For these values we can put p = 0 in Eq. (39). By making use 
formulae (36), (41), after the integration over z, the asymptotic part is presented in the form 


,^{as) 

R 


(s) = 


TTS 


27r 


■ cos ■ 


N 

E 

q=0 


(-!)»( 


TH 

a 


q—s 9 


j=0 


B Ij + 


1 — s g + s 




q + s 


(42) 


with the beta function B{x,y). In formula (42) 


Csn-I (z) == y A [(^ + d/2 - 1)2 + bD 
1=0 


(43) 


is the zeta function for a scalar field on the spacetime R x ^ and 


bu = ({D - 2){D 2)2/4 + m^r 


2^2 

H- 


(44) 


This function is well investigated in literature (see, for example, [17]) and can be presented as a 
series of incomplete zeta functions. Here we recall that the function (( 50 - 1 ( 2 ;) is a meromorphic 
function with simple poles at 2 ; = (H — l)/2 — j, where j = 0,1,2,... for D even and 0 < j < 
(D — 3)/2 for D odd. For D even one has C,gD-i{—j) = 0, j = 1,2,.... In (42), the pole term in 
the (7 = 0 summand comes from the pole of the beta function, whereas in the terms with q ^ 0 
the pole terms come from the poles of the function ^ 50 - 1 ( 2 ;). Laurent-expanding near s = 0 we 
find 

p,{as) 

Fr{s) = + F^l^ (0) + 0(s). (45) 

Using this result, for the surface energy density induced on the brane one obtains 

(46) 

where for the pole and finite contributions one has 


.(R) 

_ ^(4C 1) + 2C jp(as) 

P 

2sr^ 

R,-l’ 

.(R) 

A(4C-1) + 2C 

' p(as) 

/ 

2r^~^aSD 

-^R,0 
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+ ^^( 0 ) 


(47) 


(48) 



The corresponding formulae for the pole and finite parts of the surface stress are obtained by 
using the equation of state (21). The surface energy can be found integrating the energy density, 

^(R.surf) ^ J = ar^-^ ■ (49) 

The pole and finite parts of the vacuum expectation value of the field square on the brane are 
obtained by the formulae (22), (25), (45). 


4 Surface densities in the L-region 


In this section we consider the region between the horizon and the brane, 0 < ^ < a (L-region), 
for which one has n* = —(55^ and Kqq = —a. As in the previous section we will assume that the 
field obeys boundary condition (3) on the surface ^ = a. To deal with discrete spectrum, we 
can introduce the second brane located at ^ = 6 < a, on whose surface we impose boundary 
conditions as well. After the construction of the corresponding zeta function we take the limit 
6^0. As a result, we can see that the surface energy-momentum tensor in the L-region has 
the structure given by (19) and with the equation of state (21). For the surface energy density 
one obtains the expression 


“ = .4 = -<..4., 

2 r^ aSn 


where now Ih{A) = Fl(s)|s=o with 


1 TTS 


Fh{s) = -cos — 

'TT Z 




/=0 


Iz{ha) 


(50) 


(51) 


For a given A this expression differs from the corresponding expression for the R-region by the 
replacement Kz{x) Iz{x)- Note that the similar relation takes place for the bulk energy- 
momentum tensor as well. As in the previous section, to avoid the vacuum instability, here 
we assume that Ag < 0. Under this condition, for a given Xia the function Iz{\ia) has no real 
positive zeros with respect to z. The uniform asymptotic expansion for the integrand in (51) is 
obtained from the corresponding formula with the functions Kz{\ia) (see formula (34)) by the 
replacement 

^ (52) 

The vacuum stress is a sum of pole and finite parts 


with 


(L) = e(L) m 


+ £ 


f 


.(L) 


-(L) 


_ 4(4^ l)-|-2C ^(as) 


2sr^ 

A(4C - 1) + 2C 
2r^~^aSD 


L,-l) 


^{as) 

L,0 




( 0 ) 


(53) 

(54) 


The formulae for F^^ , Fl^^(O) are obtained from the corresponding expressions for the 

R-region by the replacements Kz{x) Iz{x) and (52). In particular. 


p(as) 

-^L 


(«) 


1 ^ / \ 
1 TTS^fruX 

-cos - > I - I 

27r 2 ^ 

9=0 


a J 


Q-s 1 / 1 _ s 

i-o ' 




CgD-l 


q + s \ 

2 J 


(55) 
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The surface energy density is related to the stress by formula (21) with the replacement R^L 
and for the total surface energy one has 

^(L,surf) (56) 


The vacuum expectation value of the field square on the brane for the L-region is also expressed 
in terms of the function For an infinitely thin brane taking the R- and L-regions together, 

the pole parts of the surface energy densities cancel for odd values of the spatial dimension D. 
In this case the total surface energy -|- £:{L,surf) -g gg^j^ j^g directly 

evaluated by the formula 


_g(surf) _ 


A{l-^C) + 2C 1 

— ' Cs-D-i ( ^ + F ) X! ^2k+i,jB ij + -,k + - 


2tt 


E 

k=0 


a J 


1=0 


+i:a 


dz 


«=o 


h{\ia) K,{Xia) ^ 2 U 2 k+i{t) 

Izika) Kz{\ia) (z^ + Afa2)fc+i 


(57) 


where = [(iV — l)/2], > 77 — 1, and t is defined by relation (41). Note that the cancellation 

of the pole terms coming from oppositely oriented faces of infinitely thin smooth boundaries takes 
place in vary many situations encountered in the literature. It is a simple consequence of the fact 
that the second fundamental forms are equal and opposite on two faces of each boundary and, 
consequently, the values of the corresponding coefficient in the heat kernel expansion summed 
over two faces of each boundary vanishes. 

We have investigated the surface densities for both R- and L-regions. In the corresponding 
braneworld scenario the geometry is made up by two slices of the region 0 < ^ < a glued together 
at the brane with a orbifold-type symmetry condition analogous to that in the Randall-Sundrum 
model (see, for instance, [6]). For an untwisted scalar field the coefficient Ag in the boundary 
condition is related to the brane mass parameter c of the field and the extrinsic curvature of 
the brane by the relation Ag = {c — C,/a)/2. For a twisted scalar Dirichlet boundary condition 
is obtained. It should be noted that in the orbifolded version due to Z 2 symmetry the extrinsic 
curvature tensor is the same on both sides of the fixed point and the cancellation of the pole 
terms for odd values D does not take place. A natural way to deal with surface divergences is 
to consider more realistic brane models with finite thickness. As it has been discussed in [18] 
for de Sitter brane model, the finite thickness of the brane regularizes the ultraviolet behavior 
and acts as a natural cutoff. 


5 Conclusion 

In this paper we have investigated the surface Casimir densities induced on a spherical brane 
in the Rindler-like spacetime Ri x by quantum fluctuations of a scalar field with an 

arbitrary curvature coupling parameter. The corresponding volume vacuum expectation values 
of the energy-momentum tensor were investigated in [6]. We consider a scalar field with Robin 
boundary conditions and as a regularization method the zeta function technique is employed. 
The spherical brane divides the background space into two regions, referred as R- and L-regions. 
We have constructed an integral representations for the corresponding zeta functions in both 
these regions, which are well suited for the analytic continuation. Subtracting and adding to the 
integrands the leading terms of the corresponding uniform asymptotic expansions, we present the 
corresponding functions as a sum of two parts. The first one is convergent at the physical point 
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and can be evaluated numerically. In the second, asymptotic part the pole contributions are given 
explicitly in terms of the zeta function for a scalar field on the spacetime R x S^~^. The latter is 
well-investigated in literature. As a consequence, the vacuum expectation values of the surface 
energy-momentum tensor for separate R- and L-regions contain pole and finite contributions. 
The remained pole term is a characteristic feature for the zeta function regularization method 
and has been found for many other cases of boundary geometries. For a minimally coupled scalar 
field, the surface energy-momentum tensor induced by quantum vacuum effects corresponds to 
a source of a cosmological constant type located on the brane. In odd spatial dimensions in 
the case of an infinitely thin brane, taking the R- and L-regions together, the pole parts of the 
surface vacuum energies cancel. As a result the total surface energy is finite and is determined 
by formula (57) with the function Uq{t) is dehned by relation (35). The results obtained here 
can be applied to the braneworld in the AdS black hole bulk in the limit when the brane is 
close to the black hole horizon. In this paper we have considered the surface energy-momentum 
tensor on a codimension one smooth brane. For non-smooth boundaries an additional part in 
the energy-momentum tensor arises located on corners. The corresponding corner terms can 
be important in codimension two braneworld scenarios (see, for instance, [19] and references 
therein). 
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